CMB Constraints on Principal Components of the Inflaton Potential 



O 
(N 



o 
u 

6 



> 

p 

O 
O 



X 



Cora Dvorkin^'^ and Wayne Hu^''^ 

' Kavli Institute for Cosmological Physics, Enrico Fermi Institute, University of Chicago, Chicago, IL 60637 

"^Department of Physics, University of Chicago, Chicago, IL 60637 

^Department of Astronomy & Astrophysics, University of Chicago, Chicago, IL 60637 

We place functional constraints on the shape of the inflaton potential from the cosmic microwave 
background through a variant of the generalized slow roll approximation that allows large ampli- 
tude, rapidly changing deviations from scale-free conditions. Employing a principal component 
decomposition of the source function G" ~ 3{V' /V)^ — 2V" /V and keeping only those measured 
to better than 10% results in 5 nearly independent Gaussian constraints that may be used to test 
any single-field inflationary model where such deviations are expected. The first component implies 
< 3% variations at the 100 Mpc scale. One component shows a 95% CL preference for deviations 
around the 300 Mpc scale at the ~ 10% level but the global significance is reduced considering the 5 
components examined. This deviation also requires a change in the cold dark matter density which 
in a flat ACDM model is disfavored by current supernova and Hubble constant data and can be 
tested with future polarization or high multipole temperature data. Its impact resembles a local 
running of the tilt from multipoles 30-800 but is only marginally consistent with a constant run- 
ning beyond this range. For this analysis, we have implemented a ~40x faster WMAP7 likelihood 
method which we have made publicly available. 



I. INTRODUCTION 

Under the assumption that cosmological perturbations 
were generated during an inflationary period from quan- 
tum fluctuations in a single scalar field, features in the 
cosmic microwave background (CMB) temperature and 
polarization power spectra constrain features in the in- 
flaton potential V {(/)). 

The usual slow-roll assumption of a small and nearly 
constant dimensionless slope {V /V)'^ and curvature 
{V" /V) of the inflaton potential leads to featureless 
power law initial curvature fluctuations. Breaking any 
of these assumptions leads to features in the initial spec- 
trum. Indeed glitches in the observed temperature power 
spectrum of the CMB [l| have led to recent interest in 
exploring models with strong features in the potential 
{e.g. [2i9|). The significance at which such models are 
favored is difficult to address due to the a posteriori man- 
ner in which they were chosen. 

Here we take a more model independent approach to 
constraining the shape of the inflaton potential. We have 
recently shown that even in the presence of large local 
changes in the curvature of the inflaton potential that 
can explain the glitches, there is to excellent approxima- 
tion only a single function of the inflaton potential that 
the observations constrain Q- Moreover, this function 
is approximately the same combination of slope and cur- 
vature that enters into the calculation of the scalar tilt 
in the ordinary slow-roll approximation. In this gener- 
alized slow-roll (GSR) formalism this quantity need not 
be small or constant |10l - [l3 |. With it, one can bypass 
a parameterization of the initial curvature power spec- 
trum {e.g. |13l - [l9| ) and the problem that not all spectra 
correspond to possible inflationary models. 

In this Paper, we take a principal components ap- 
proach |20| to functional constraints on the inflaton po- 
tential under the GSR formalism. Principal compo- 



nents constructed a priori from a noise model of the 
WMAP CMB measurements determine the theoretically 
best constrained deviations from a featureless potential 
before examining the actual data. Constraints from the 
low order principal components thus efficiently encap- 
sulate the expected information content of the data and 
may be used to test a variety of inflationary models with- 
out a reanalysis of the data. 

We begin in SJIT] with a brief review of the GSR for- 
malism. In ijIIIBI we develop the principal component 
analysis of the GSR source function and apply it to the 
WMAP 7 year (WMAP7) data in f|TVl In ^we consider 
applications of these derived constraints on the inflaton 
potential and discuss these results in §VII In the Ap- 
pendix, we present the fast likelihood approach to the 
WMAP7 data employed in these analyses. 



II. GENERALIZED SLOW ROLL 

Given a specific model for the inflaton potential, the 
initial curvature fluctuation spectrum can always be nu- 
merically computed by solving the linearized scalar fleld 
equation. The slow roll approaches, however, provide 
model independent mappings from the inflaton potential 
to the curvature spectrum provided that the requisite 
approximations hold. 

The generalized slow roll (GSR) formalism was orig- 
inally developed to consider cases where the usual slow 
roll parameters en = {(f>/H)'^ /2 and tjh = —ilHcj) are 
small but rjn is not necessarily constant [l(]|4l2| . Here 
<j) is the inflaton field and overdots represent derivatives 
with respect to coordinate time t. In a previous paper 
3 , we have shown that a variant of GSR works well for 
cases where tjh becomes of order unity for a fraction of 
an e-fold. 

In this variant of the GSR approximation, the curva- 



ture power spectrum is a non-linear functional 
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which is related to the inflaton potential through the 
background solution / = 2TT(j)ari/H. Primes here and 
below denote derivatives with respect to In 77 and rj = 
/j ""'^ dt' /a is the conformal time to the end of inflation. 
We require fcmaxf?min < 1 and ki^inVmax > 1 for the range 
in k that we are calculating the spectrum. Integrating G" 
gives 
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define the linear and nonlinear response of the curvature 
spectrum to G' respectively. For the models we consider, 
the nonlinear response is small compared with the linear 
one. 

The key property of Eq. ([1]) is that deviations from 
scale invariance in the curvature spectrum depend only 
on a single function of time G". Moreover, to good ap- 
proximation, this function is related to the inflaton po- 
tential as M 



G'«3(^)2-2^ 
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so long as jry^j ^ \vh\ when tjh is large, i.e. that rj^ 
remains large only for a fraction of an e-fold Q . Finally, 
if the ordinary slow roll approximation where en and 
r]H are both small and nearly constant holds, then G' 
may be evaluated at horizon crossing 77 « k~^ and taken 
out of the integrals in Eq. ([T]). As we shall see below, 
under this approximation G' = 1 — ris. By allowing G' 
to be both time varying and potentially large, we recover 
ordinary slow roll results where they apply but allow the 
data themselves to test their validity. 



III. PRINCIPAL COMPONENTS 

The GSR approximation allows us to go beyond spe- 
ciflc models of inflation in examining how the data con- 
strain the inflaton potential. The data directly constrain 



the function G' and hence the derivatives of the inflaton 
potential through Eq. ([5]). 

Given that G' is related to the curvature spectrum A^ 
by an integral relation and the curvature spectrum itself 
is related to the observable CMB power spectra by a line- 
of-sight integration, not all aspects of the function G' are 
observable even with perfect data. 

We therefore seek a basis for an efficient representation 
of observable properties of G'. We begin in ^III Al with 
a general description of a basis expansion for G' and its 
relationship to the usual normalization and tilt param- 
eters. We then turn in ijIIIBI to principal components 
(PCs) as the basis which best encapsulates expected de- 
viations from scale-free conditions 1201 • 



A. Basis Expansion 

In general, we seek to represent the function G' as 

G'(lnr/) = ^m„5,(ln77), (6) 



a=0 

where the basis functions Sa for a > are assumed to 
have compact support in some region between In 77niin and 
ln'7max corresponding to the range probed by the data. 
We assume S'o = 1 so that within this range niQ repre- 
sents a constant tilt in the spectrum. 

We seek functions that obey the orthogonality and 
completeness relations 
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where A In 77 = Inr^max — Inrj^ii^. From these relations, 
the TTia amplitudes are related to G' as 
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Note that our normalization differs from that of Ref. [20| 
in that unit amplitude 777,0 corresponds to unit variance 
in G' averaged over the whole range in In 77. 

Substituting this model into the power spectrum ex- 
pression ([1]) yields 



lnA|.(fc) w G(ln77,„in) + 777oG 
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l + ^[E'"aXa(fc)j 

where G = |-7E-ln2« 1.06297 with 7e as the Euler- 
Mascheroni constant, Xq — 7r/2. The fc-space responses 
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FIG. 1: The first 5 principal components (PCs) of G' as a function 
of conformal time based on the WMAP7 specifications. The power 
law model with zero amplitude PCs is shown in red dashed lines. 
The first 5 PCs represent a local expansion of G' around -q ~ 10^ 
Mpc. 



FIG. 2: Predicted RMS error on the PC amplitudes as a function of 
mode number for WMAP7 data. PCs are rank ordered from lowest 
to highest error with the first 5 describing modes with better than 
~ 10% constraints on G' . 



to the modes are characterized by 
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where /cTymin ^ 1 and krjYanx ^ 1- Note that if jtiq = 
for a > 

from which we infer that the model is a pure power law 
spectrum. We can therefore choose instead to represent 
G'(ln?7,„in) and ttiq by Us and As bringing our parameter- 
ization of the power spectrum to 
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Note that this replacement ensures that the normaliza- 
tion and tilt parameters are defined at a scale kp that is 
well-constrained by the data. Hence any small and un- 
observable running of G and G" from ?7inin to ry ~ 1/kp 
is absorbed into As and n^ respectively. Non zero values 
of 7tJq>o required by the data thus represent a deviation 
from purely scale-free initial conditions. 

We hereafter consider Eq. (fTTj) as the definition of the 
parameterized curvature power spectrum. In practice we 
choose kp = 0.05 Mpc~^, and for reference note that k sa 
0.02 Mpc~^ for modes contributing to the well-measured 
first acoustic peak. 



B. Principal Component Basis 

We choose here to construct the basis functions Sa 
from the principal components (PCs) of the projected 
WMAP7 covariance matrix for perturbations in G". To 
define the PCs, we begin with a fiducial fiat ACDM model 
with a scale-free initial spectrum. We take the baryon 
density to be i^^h'^ = 0.02268, cold dark matter den- 
sity flch^ = 0.1080; cosmological constant Ha = 0.7507, 
optical depth t = 0.089, As = 2.41 x 10-^ Us = 0.96 
and rriayo = 0. This model corresponds to a constant 
G' = 0.04. 

We then construct the PCs as the theoretically best 
constrained non-constant deviations in G" around this 
fiducial model. We start by adding a set of perturbations 
in G'(ln77j) = 0.04 + pi at 50 equally spaced intervals in 
In 77 between 1 < ?7/Mpc < 10^. From this set, we define 
the function G'(ln77) by a cubic spline. This sampling of 
10 per decade or Slnrj = 0.23 across A In 77 = 11.5 is suf- 
ficient to capture the observable properties of G' barring 
unphysical models with both high frequency and high 
amplitude perturbations. The spline ensures a smooth 
interpolation between the samples. 

To these parameters pi we add the cosmological pa- 
rameters 



Pm = {pi,As,ns,T,nbh'^,flch'^,9} . 
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where the angular extent of the sound horizon 9 takes 
the place of Q\, and construct the Fisher matrix 
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FIG. 3: Sensitivity of the curvature power spectrum to the first 5 
PC parameters. Low order PCs mainly change the power spectrum 
at wavenumber in the decade around k ~ 10~^ Mpc~^. Red dashed 
hne represents the zero PC amphtude fiducial power law model. 



FIG. 5: Predicted RMS errors on running of tilt a as a function 
of the maximum number of PC components included. Note that 
the errors cease to improve after the first 4 components and most 
of the improvement comes from the 4th component. 
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FIG. 4: Sensitivity of the temperature power spectrum to the first 
5 PC parameters. Low order PCs represent slowly varying features 
around the first peak i ~ 200. Red dashed line represents the zero 
PC amplitude fiducial power law model. 
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where the WMAP7 noise power N^ 
inferred from the temperature power spectrum errors 
from the LAMBDA site [23| and the assumption that 
Nf^ ^ 2Nj^. We set f^ax = 1200 and /sky = 1. 

We then invert the Fisher matrix to form the covari- 



ance matrix C 



Next we take the sub-block Ci 



and decompose it with the orthonormal matrix Sja, 
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rank ordered from lowest to highest a a ■ Each eigenvector 



defines a discrete sampling of the basis function Sa via 
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with the normalization of Eq. ([7|). The full functions 
Sa (In 7/) arc defined by taking a cubic spline through the 
samples. The first 5 PC components are shown in Fig. [1] 
In the Fisher approximation aa represents the 
WMAP7 expected errors for uia 



{marrib) = 5ab<yl , 



(17) 



for a zero mean fiducial model {ma) = 0. 

These errors arc shown in Fig. [2l Note that the noise 
rises by a factor of 4 across the first 5 components and 
then increases to order unity by the 20th component. 
Given that the peak amplitudes of Sa lie in the ~ 4 range, 
a ^ 2.5% error corresponds to a ^ 10% peak variation 
in G' or the effective tilt. By keeping the first 5 PCs, we 
retain all of the constraints on deviations in G' in the 10% 
range. In this sense, these 5 modes represent the best 
set for examining deviations from a scale-free power law 
model. Note that no actual WMAP7 data goes into the 
construction and so that the modes are chosen a priori 
rather than a posteriori. 

For the first 5 PC components, the Sa basis functions 
arc centered near 7/ w 10^ Mpc and rcfiect the strong 
WMAP sensitivity to the first peak at ^ sa 200 or A: - 0.02 
Mpc~^ in the fiducial model. The first 5 components 
resemble a local decomposition of G" in the decade sur- 
rounding this scale. This fact can be seen more directly 
by examining the sensitivity of the curvature and tem- 
perature power spectra to the 5 PC amplitudes (see Figs. 
[3]and|31 respectively). 

As an example of the utility of retaining only the first 
5 PCs, consider a linear model for G" 
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The local slope in the power spectrum in the GSR ap- 
proximation is 



dlnA^ 
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and so a = dus/dlnk and is equivalent to the running 
of the tilt. This linear model can be projected onto the 
first 5 PCs. The variance of a is then given by 
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By virtue of the marginalization of n^ in the construc- 
tion of the PCs, the Sa functions have nearly zero mean 
and dma I da does not depend on the scale rjQ where the 
effective tilt is defined as G'(ln7yo) = 1 — n^. 

Fig. [S] shows the predicted RMS error on a as a func- 
tion of the maximum PC mode retained. With 5 PCs, 
a{a) = 0.0355, while the fully-saturated value of the er- 
ror with 50 PCs is (j{a) = 0.0327. This should be com- 
pared with the projected error using a itself as a Fisher 
matrix parameter, (T(a) = 0.0328. These results verify 
the completeness of the 50 PC basis as well as show that 
most of the information for |a| <C 1 is expected to come 
from the first 5 PCs. 

In fact most of the information comes from a single 
mode, the 4th. This mode corresponds to a local varia- 
tion in the effective tilt G" with a null near 300Mpc (see 
Fig. [ij, or in the power spectrum near ^ 0.003 Mpc~^ 
(see Fig. [3]), and an extent spanning 1-2 decades. Even 
though the 1st mode has smaller overall errors and bet- 
ter constrains peak variations in G", it is not the most 
effective mode for constraining running of the tilt given 
its extremely local form. 

This caveat applies more generally. A given model may 
have large deviations in G' in a region where the data 
does not best constrain G". In this case the first 5 PCs 
no longer represent a complete or efficient representation. 
We return to this point in [JV] 



C. MCMC 

We use a Markov Chain Monte Carlo (MCMC) likeli- 
hood analysis to determine joint constraints on the first 
5 PC amplitudes and cosmological parameters 



Pa* == {"^1 
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,T,ni,h\n,h\0}. (22) 



On top of this basic set we also examine the impact of 
spatial curvature Hk and tensor-scalar ratio r and the 
amplitude of a Sunyaev-Zel'dovich contaminant Asz on 
a case-by-case basis. 



The MCMC algorithm samples the parameter space 
evaluating the likelihood jC{x\p) of the data x given each 
proposed parameter set p {e.g. see [2l|,[22|). The poste- 
rior distribution is obtained using Bayes' Theorem, 



7'(p|x) 
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where V{p) is the prior probability density. We place 
non-informative tophat priors on all parameters in 
Eq. (|22|) . For example, for the PC amplitudes we take 
'P{ma>Q) — 1 for —1 < ma>o ^ 1 and otherwise. 

The MCMC algorithm generates random draws from 
the posterior distribution that are fair samples of the like- 
lihood surface. We test convergence of the samples to a 
stationary distribution that approximates the joint pos- 
terior density 7^(p|x) by applying a conservative Gelnian- 
Rubin criterion [2J] ofi?— 1<0.01 across four chains. 
We use the code CosmoMC ^ for the MCMC analysis 

m- ^ 

For the WMAP7 data [25[, we optimize the likelihood 
code available at the LAMBDA web site as detailed in 
the Appendix. The net improvement in speed on an 8- 
core desktop processor is a factor of ^ 40 which should 
enable future studies with 20 or more PCs or other initial 
power spectrum parameters. 



IV. CONSTRAINTS 

In this section we present the constraints on the 5 best 
measured principal components of the GSR source func- 
tion G" and implicitly the inflaton potential V {(/)). We 
first examine constraints using the WMAP7 data only 
and then joint with a variety of cosmological constraints 
to remove residual parameter degeneracies. 



A. WMAP 

We begin by considering the WMAP7 data in a fiat 
ACDM cosmological context. We show the probability 
distributions of the first 5 PCs and cosmological param- 
eters in Fig. [5] Table H] shows the mean, standard devia- 
tion of the posterior probabilities as well as the maximum 
likelihood parameter values for the power law models vs. 
the 5PC models. 

For visualization purposes we show in Fig. [7] the func- 
tional posterior probability of 



G'^{ri) = 1 - rts + ^ maSa{r]) 



(24) 



which should be interpreted as G" filtered through the 
first 5 PCs and not a reconstruction of G" itself ( cf. Fig.fTT] 
below). 

All 5 PCs are tightly constrained, with errors that 
are comparable to the Fisher matrix projection, nearly 
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FIG. 6: Posterior probability distributions of the cosmologieal and 5 PC parameters using WMAP7 data. The red dashed hne represents 
the power law results where the first 5 PC parameters are held fixed to ma = 0. 
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FIG. 7: The 5 PC filtered G'^ posterior using WMAP7 data. The 
shaded region encloses the 68% CL region and the upper and lower 
curves show the upper and lower 95% CL limits. The maximum 
likelihood (ML) Gg is shown as the thick central curve, and the 
power law ML model is shown in red dashed lines. Structure in 
this representation mainly reflects the form of the PC modes and 
is dominated by the modes with the largest uncertainties. 



Gaussian posteriors and little covariance with each other. 
The correlation coefficients between two different iria's, 
|C'r„^m[,/c(?nQ)cr(7Ti(,)| < 0.2. The first component mi in 
particular is consistent with zero and places the tight- 
est constraints of < 3% local variations in G" around 
77 w 10^ Mpc. Interestingly, the power law prediction of 
7714 — lies in the tails of the posterior with as extreme 
or more values disfavored at 94.8%CL. With the Gaus- 
sian approximation 7714 — is 1.9(t from the mean and 
disfavored at 94.5%CL. 

In Fig. [HI we show joint posteriors of the PCs with 
other parameters. Notably, for the anomalous 777,4 com- 
ponent there is a degeneracy with Q.ch?', flbh'^ and 77s 
which is also refiected in the broadening of the cosmo- 
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FIG. 8: Joint probability distributions of the principal component 
amplitudes and the cosmologieal parameters from MCMC analysis 
of WMAP7 data (68% and 95% CL contours). 



logical parameter posteriors in Fig. |S] and the shift in 
means and maximum likelihood values in Table HI The 
maximum likelihood (ML) model found by the chain is 
an improvement over the power law case of 2 A In £ w 5 
which is marginal considering the addition of 5 parame- 
ters. 

The intriguing aspect of the ML model, like the uia 
posteriors, is that the improvement is concentrated in 
the TO4 component. Note that the finite 7774 component 
allows 77s = 1 to be a good fit to the data implying that 
the data can be marginally better fit by a local deviation 
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FIG. 10: Decomposition of the fractional difference between the 
PC and power law (PL) maximum hkehhood models (ML) shown 
in Fig. |9] into contributions from specific parameters. Curves show 
the cumulative effect of adjusting the PL ML parameters to the PC 
ML values. The main effects come from the change in nn, Q.ch? , 
and Us. 
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ization cross correlation part of the likelihood (MASTER 
TETE), where 2A\nC^ 3.5, with half of this contribu- 
tion coming from £ < 200. The low-.^ temperature part 
of the likelihood has an improvement of 2Aln£ =1.8 
relative to the PL model, and there is a smaller improve- 
ment coming from the high-£ (MASTER TT) part with 
2Aln£ =1.1. Finally, the PC model is worse than the 
PL model by 2 A \nC = —1.3 for the low i < 24 polariza- 
tion. 



FIG. 9: Power spectra of the 5 PCs maximum likelihood model 
(black curve) compared to power law (PL) maximum likelihood 
model (red dashed curve) in the upper panel, and the difference 
(PC-PL)/PL in the lower panel. Top: temperature power spec- 
trum. Bottom: polarization power spectrum. 



from scale invariance rather than tilt. 

An examination of the ML model helps illuminate 
the degeneracies with cosmological parameters. Fig. [9] 
(top) shows the temperature power spectra of the 5 PC 
ML model compared to the power law ML model (up- 
per panel), and the fractional difference between the 
two (lower panel). Note that in the well constrained 
£ ~ 30 — 800 regime the two spectra agree to ~ 1% 
or better. This accidental degeneracy is not preserved 
beyond i = 1000. Furthermore the E-mode polarization 
power spectra shown in Fig. [5] (bottom) reveal substan- 
tially larger fractional deviations of up to ^ 10% that 
break the degeneracy in the temperature spectra. 

Indeed the main improvement of the PC model relative 
to the PL model actually comes from the £ > 2A polar- 



The intermediate £ temperature degeneracy exhibited 
by the two models is further explored in Fig. [TOl We 
show the impact of the cumulative parameter variations 
between the PL and PC ML models. The addition of 7714 
carries almost all of the impact of the PCs and is mainly 
compensated by variations in ftch^ to adjust the relative 
heights of the peaks and n^ to tilt the spectrum to small 
scales. The change in ilch^ also changes the physical size 
of the sound horizon which must be compensated by a 
change in the distance to recombination reflected in a 
lower value for Hq and Ha to leave the angular scale 9 
compatible with the data. 

These results are robust to marginalizing Qk given any 
reasonable prior on Hq, an SZ component through Asz, 
or tensors within the _B-mode measured limits of the BI- 
CEP experiment r = 0.03to:26 111- Through an explicit 
MCMC analysis of these separate cases, we have veri- 
fied that the shift in the means and change in the errors 
for all 5 PCs are much smaller than la. The largest 
effect is from marginalizing tensors where for example 
TO4 = 0.042 ± 0.20. A small improvement in the i?-mode 
limits would eliminate this ambiguity entirely. 
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2.217 


2.040 ±0.196 


2.067 


n,h' 


0.1116 ±0.0053 


0.1130 


0.1308 ±0.0127 


0.1284 


e 


1.0386 ± 0.0026 


1.0387 


1.0361 ± 0.0049 


1.0365 


T 


0.088 ±0.014 


0.088 


0.089 ±0.017 


0.088 


Us 


0.9650 ±0.0136 


0.9622 


0.9916 ±0.0233 


0.9877 


lnlW'"As] 


3.083 ± 0.034 


3.088 


3.119 ±0.041 


3.105 


mi 








0.0014 ± 0.0077 


0.0021 


7712 








0.0015 ±0.0132 


0.0068 


7713 








-0.0253 ±0.0197 


-0.0264 


7714 








0.0339 ±0.0175 


0.0337 


77l5 

Ho 








-0.0033 ±0.0315 


0.0023 


70.13 ± 2.38 


69.50 


61.41 ±6.08 


62.18 


Oa 


0.726 ± 0.028 


0.720 


0.581 ±0.116 


0.614 


-21n£ 


7474.97 


7469.82 



TABLE L Means, standard deviations (left subdivision of columns) and maximum likelihood values (right subdivision of columns) with 
likelihood values for ACDM and the 5 PCs model with WMAP7 data. 



Joint Constraints 



APPLICATIONS 



The results of the previous section suggest that other 
data which measure the high-£ temperature spectrum, 
polarization spectrum, or pin down the cosmological pa- 
rameters that control the distance to recombination and 
baryon density can eliminate the remaining degeneracies 
and enable WMAP7 to better constrain the inflaton po- 
tential. 

We start with adding more CMB information from the 
QUAD experiment. QUAD helps mainly by reducing the 
7715 ~ ^bh'^ degeneracy. Interestingly, most of the impact 
comes from the polarization measurements rather than 
the extended range of the temperature constraints. 

Adding in non-CMB cosmological information helps 
even more, especially with 7^4. We take the UNI0N2 
supernovae dataset |3l|, the SHOES Hq = (74.2 ± 3.6) 
km/s/Mpc measurement |27| and a big bang nucleosyn- 



thesis constraint of rj^/i^ = 0.022 ± 0.002. In a flat 
ACDM cosmology, the degeneracy between 7774 and flch^ 
is nearly eliminated yielding 7774 = 0.0191 ± 0.0163, 
i.e. consistent with power law initial conditions. Of the 
additional data it is the supernovae that drive this im- 
provement by disfavoring the low J7a values required by 
the increase in il.ch^. 

On the other hand, these improvements require an as- 
sumption that the dark energy is a cosmological constant 
and the Universe is spatially flat. For example if flu is 
marginalized, 7774 = 0.0384±0.0197. The addition of spa- 
tial curvature allows the freedom to adjust the relative 
distance to the high-z supernova and recombination. A 
better measurement of Hq could resolve this degeneracy 
since the constraints on fix are already limited by the 
SHOES data. 

Table [ll] summarizes these results for the constraints 
on the PC amplitudes. 



The model independent results of the previous section 
can be used to test a wide variety of inflationary devia- 
tions from scale-free conditions. Moreover given that the 
constraints on the PC amplitudes are uncorrelated and 
approximately Gaussian, these tests are straightforward 
to apply. 

As the simplest example, take the running of the tilt 
model defined by a = dris/dlnk through Eq. (fT5)) . Using 
Eq. ([8]), we obtain 

7771= 0.048a, 7772 = -0.079a, 7773 = 0.054a, 
7714 = -0.576a , 7775 = -0.034a . (25) 

We can then construct the effective x^ statistic 



x^(«) = E 



77la(a) 



(26) 



With the means and variances taken from Table |T] for 
WMAP7 we obtain a = -0.057 ± 0.029. When we take 
into account the covariance between the PC amplitudes, 
we obtain a 3% shift in the mean with the same error: 
a = —0.059 ±0.029. Likewise we have verified that using 



Eq. (p5|) in a separate MCMC analysis gives consistent 
results a = -0.058 ± 0.030. 

This result should be compared with the direct analysis 
of the running of the tilt which gives a = —0.034 ± 0.027 
consistent with the analysis from [25|. The mean is 
shifted from the PC derived mean by ^ O.Sa while the 
errors are 7% higher. 

The small shift in the mean is driven mainly by the 
truncation to 5 PC components. This is in spite of the 
Fisher expectation in i jHIBI that for an infinitesimal a 
the first 5 PC components contain nearly all the infor- 
mation. A running of the tilt model where a = —0.057, 



PC 


+QUAD 


+BBN+SN+iyo, flat 


±BBN±SN±//o, w/Qa- 


mi 


0.0000 ± 0.0072 


0.0045 ± 0.0071 


0.0027 ± 0.0073 


m2 


0.0033 ± 0.0123 


0.0091 ±0.0121 


0.0045 ± 0.0125 


m-i 


-0.0261 ± 0.0184 


-0.0120 ± 0.0166 


-0.0208 ± 0.0178 


7714 


0.0296 ± 0.0168 


0.0191 ±0.0163 


0.0384 ± 0.0197 


m,5 


0.0149 ± 0.0250 


0.0187 ±0.0249 


0.0091 ± 0.0256 



TABLE II: Means and standard deviations of the posterior probabilities of the PC amplitudes with different data sets added to the 
WMAP7 data. For supernovae (SN) we used the UNION2 dataset, and for Ho the SHOES measurement. 



which fits the intermediate ^ ^ 30 — 800 range well, im- 
plies a large change across the extended observable range 
from ?7 - 20 - 5000 Mpc of \5G'\ - |<5n,| - 0.3. In 
particular, it overpredicts the suppression of the ^ < 30 
temperature multipoles. For the same amplitude and tilt 
at the first peak, the amplitude at the horizon is sup- 
pressed by ~ e^"/^^'" ^^°°) = 0.55. This suppression can 
only be partially compensated by red tilting the spec- 
trum without over suppressing the high (. > 800 multi- 
poles. Note that cosmic variance completely dominates 
the uncertainties in the £ < 30 region and decreases with 
the predicted signal, an effect that is not represented in 
the Fisher matrix. 

In other words, the 2cr preference for finite 777.4 is not 
completely consistent with a constant running of the tilt 
but rather points to a more local deviation from scale- 
free conditions. When we eliminate this preference by 
adding in the additional SN, Hq and BBN constraints 
the inferred limits on a from the first 5 PC amplitudes 
becomes a = —0.033 ± 0.027. Note however that the 
direct constraints on a also improve with the addition of 
these data sets to a = -0.013 ± 0.021. 

In Fig. [TT] we plot an a = —0.033 model for the 5 
PCs filtered 6*5 against the posterior constraints from 
the WMAP7 data and additional SN, H^ and BBN con- 
straints in a flat universe. We overplot the original unfil- 
tered G" for this a and note that even with the reduced 
value the deviations become large outside of the region 
probed by the first 5 PCs. 

This example points out a caveat to the use of the first 
5 PCs as general constraints on models. For a model 
with features that are substantially larger in a regime 
away from the well-constrained first acoustic peak, the 
first 5 PCs may not be the best constraints in terms of 
signal-to-noise. One can check whether this is the case by 
examining the predicted G" or by projecting the model 
onto the full 50 PC space and checking for large ampli- 
tude components. Indeed if the higher components are 
extremely large compared with the low components, non- 
linear effects can break the orthogonality of PCs and lead 
to larger allowed variations in the low components when 
compensated by the high components. 

As an example, the full 50 PC decomposition of the 
the step function potential model from Q which fits the 
£ = 20 — 40 glitches in the temperature power spectrum is 
shown in Fig. [T^l Interestingly, 7774 = 0.0266 in the step 




77 [MpcJ 



FIG. 11; The 5 PC filtered G'^ posterior using WMAP7 data 
and additional SN, Hq and BBN constraints in a flat universe. 
The shaded region encloses the 68% CL region and the upper and 
lower curves show the upper and lower 95% CL limits. A model 
with running of the tilt a = —0.033, the mean value given these 
constraints, is shown as the thick solid blue curve and the ML PL 
model as the dashed red curve. The unfiltered G' of the same 
a is shown in blue dashed lines for comparison (arbitrary offset). 
Note that outside the range probed by the first 5 PCs the model 
deviations continue to grow linearly and oscillating features in G'r 
do not necessarily imply features in the underlying G' . 



model and has the highest amplitude of the first 5 com- 
ponents. On the other hand, a complete analysis based 
on signal-to-noise would require ~ 20 PC components. 
By keeping only 5 components, the improvement com- 
pared with the ML PL model is only A^-^ = —1.7. In 
other words, the step function model is certainly allowed 
by our 5 PC constraint and even marginally favored but 
the majority of the improvement is not captured by the 
truncated analysis. 

Nevertheless, when interpreted as an upper limit on 
deviations from scale-free conditions, the 5 PC approach 
works as a general, albeit typically conservative, method 
to constrain a wide variety of possible deviations from a 
single analysis. As the running of the tilt example shows, 
the results are remarkably close to a direct analysis and 
differences can be used to expose the self-consistency of 
the model inferences with independent parts of the data. 
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FIG. 12: Principal component amplitudes for the step function 
potential model [5|| that best fits the glitches in the temperature 
spectrum at £ ~ 20 — 40 (upper panel), and projected cumulative 
{S/N)'^ = ^(ma/ua)^ (lower panel). Given the large values of 
ma in the high order PC components, ~ 20 PCs are required to 
completely characterize this model. 



VI. DISCUSSION 

We have employed a variant of the generahzed slow roll 
approximation (GSR) that allows large amplitude devi- 
ations from scale-free conditions to place functional con- 
straints on the GSR source function and implicitly on the 
inflaton potential. By employing a principal component 
(PC) decomposition, we isolated the 5 best functional 
constraints imposed by the WMAP7 data. The analysis 
is greatly facilitated by our optimization of the WMAP7 
likelihood code which we have made publicly available 

These 5 PCs provide incisive constraints on the infla- 
ton potential around the e-folds of inflation when the 
scales associated with the first acoustic peak were cross- 
ing the horizon. Non-zero values for their amplitudes 
represent deviations from slow roll and power law initial 
spectra. The first component implies that deviations are 
less than ~ 3% near ij ~ 10^ Mpc and the first 5 rep- 
resent constraints around that scale at better than the 
10% level. The result is 5 nearly independent Gaussian 
constraints that can be applied to any inflationary model 
where this level of deviation is expected. We have also 
made the eigenfunctions, which are required to project a 
given model onto the PC amplitudes, publicly available. 
These limits are robust to the inclusion of tensor contri- 
butions allowed by current B-mode limits, spatial cur- 
vature and Sunyaev-Zel'dovich contamination from un- 
resolved clusters. 

Interestingly, for the 4th principal component the null 
prediction of scale-free initial conditions is disfavored at 
the 95% CL. However, given the 5 added parameters, this 



result does not rule out a power law initial spectrum at 
a signiflcant level. Moreover, the relatively large devia- 
tions implied by this anomalous mode are allowed only 
through correspondingly large variations in the cosmo- 
logical parameters, mainly the cold dark matter and its 
effect on the sound horizon and by inference the distance 
to recombination. 

Further information from the CMB polarization and 
high-£ temperature power spectrum can break this de- 
generacy. The QUAD polarization data already have 
some impact on the constraint and the Planck satellite 
should deflnitively resolve this issue. External data also 
can break this degeneracy. In particular in a flat ACDM 
cosmology, the distance to high redshift supernovae re- 
duce the preference for finite m^ from 1.9cr to 1.2cr. How- 
ever this improvement disappears if spatial curvature is 
marginalized. 

This anomalous 4th PC resembles a local running of 
the tilt around scales of j] ^ SOOMpc. Direct analysis of 
a global constant running of the tilt shows that this pref- 
erence is mainly local, i.e. the low and high multipoles 
prefer a different and smaller running than the intermedi- 
ate multipoles that the first 5 PCs probe. The running of 
the tilt example illustrates the use of the PC constraints 
both as a technique to constrain inflationary parameters 
arising from different models with a general analysis and 
as a method for examining what aspects of the data drive 
the constraints. 

The running of the tilt example also illustrates that 
for models where deviations from scale-free conditions 
become much larger than ^ 10% away from the well- 
constrained region of the acoustic peaks, more principal 
components are required to ensure a complete and inci- 
sive description. We intend to examine these issues in a 
future work. 
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Appendix A: Fast WMAP Likelihood Evaluation 

In this Appendix, we describe the optimization of 
the WMAP likelihood code and fast approximate tech- 
niques for describing the low-i polarization information. 
Changes in the initial power spectrum do not require re- 
computation of the radiation transfer function and are so- 
called fast parameters for CosmoMC. Hence the WMAP7 
likelihood computation is the main bottleneck for the 
MCMC analysis. 

We first OpenMP parallelize the likelihood code and 
remove bottlenecks in the computation of the temper- 
ature and high-£ polarization likelihood. We obtain a 
^ 2.&Ncoj:e speedup of those parts of the likelihood where 
A^corc is the number of cores in a shared memory machine. 
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FIG. 13: Comparison of the low-^ polarization pixel likelihood 
— 21n£^j and the approximate fit as a function of _E-mode po- 
larization amplitude in two multipole bands p^i {I = i — 6), Pe'z 
(f = 8). Models from a power law chain (red crosses) and from a 5 
PCs chain (black crosses) whose likelihood relative to the minimum 
1645.84 are (within ±0.1) of 2.29 (68.27% CL) and 6.18 (95.45% 
CL) are shown with the contours from the fit overplotted (blue 
curves) . 



These changes exactly preserve the accuracy of the hke- 
hhood evaluation. 

In place of the computationally expensive low-£ po- 
larization pixel likelihood, we seek a fast but accurate 
approximation. The WMAP team has shown that most 
of the information in the power law ACDM parameter 
space lies in multipoles £ ~ 2 — 7 as essentially an overall 
amplitude of power [25|. However, in the broader param- 
eter space allowed by the PCs of G', we find that a single 
amplitude is insufficient to describe the information con- 
tent of the pixel likelihood. 

Instead we fit the likelihood function to a two band 
approximation 



PEi 



E 



2ti 



EE 



where the first band i = 1 has ^imin 
and the second band z = 2 has £2min 



1/2 



= 4, h. 



2 max 



£■ 



(Al) 



= 6 
and 



Ae, 



- 1 normalizes the parameter to reflect 



the average bandpower. We find that the pixel likelihood 
is well approximated by a Gaussian in these two bands 
for the models under consideration 

- 21n4<24 ^A+ipE- PEfC-\pE - Pe) (A2) 

with the parameters pEi = 0.2614 fiK, pE2 — 0.01955 




-0.04 -0.02 0.00 0.02 0.04 0.06 0.08 0.10 



FIG. 14: Posterior probability distribution of the optical depth 
T and the fourth PC amplitude using the exact likelihood (black 
curves) and the approximation (blue dashed curves) with WMAP 
data only. 



^K, A = 1645.84 



498.31 -214.23 
-214.23 190.23 



fiK' 



(A3) 



In Fig. [13] we show the accuracy of the fit compared with 
the pixel likelihood for both power law models and mod- 
els with additional 5 PCs of G". Note that the power law 
models lie on a ID curve in this space and can be well 
parameterized by a single amplitude whereas the 5PC 
models do not. In fact, in the 2 band space models with 
low Pei and high pE2 that populate a direction nearly 
orthogonal to the power law models are more strongly 
constrained than the total power at low £ would suggest. 
Since this approximation has trivial computational cost, 
the net improvement in speed is approximately ^ bNcorc- 
For the cases of interest, the approximation works re- 
markably well. As an example, we have run an MCMC 
chain with the exact pixel likelihood for the 5PC chain 
with WMAP data only. In Fig. [T31 we compare the pos- 
terior probability distribution of the optical depth r and 
TO4 using the full likelihood and the approximation. The 
difference between the pixel likelihood and the approx- 
imation for the 5 PCs maximum likelihood model with 
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WMAP data is likewise negligible: |2Aln£| = 0.05. 

We have also checked that the likelihood approxima- 
tion remains valid to 10% or better in the (2Aln£)^/^ 
significance of differences between models with varying 
reionization historyas parameterized by ionization prin- 
cipal components [28J. 

Larger differences can occur for models with sharp, or- 
der unity, initial power spectrum features at the horizon 
scale. These project onto the temperature and polariza- 



tion spectra differently and leads to qualitatively different 
results for the temperature-polarization cross spectrum. 
In this case one can get discrepancies of order unity in 
2Aln£ that err on the side of allowing discrepant mod- 
els. Even these differences can typically be taken into 
account via importance sampling at a much smaller com- 
putational cost than evaluating the exact pixel likelihood 
during the MCMC run itself. 
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